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Abstract. Radial and vertical profiles are determined for a sample of 34 edge-on disk galaxies in the HDFs, selected 
for their apparent diameter larger than 1.3"and their unperturbed morphology. The thickness and flatness of their 
galactic disks are determined and discussed with regard to evolution with redshift. We find that sub-L* spiral 
galaxies with z ~ 1 have a relative thickness or flatness (characterized by h z /h the scaleheight to scalelength 
ratio) globally similar to those in the local Universe. A slight trend is however apparent, with the h z /h flatness 
ratio larger by a factor of ~ 1.5 in distant galaxies if compared to local samples. In absolute value, the disks 
are smaller than in present-day galaxies. About half of the z ~ 1 spiral disks show a non-exponential surface 
brightness distribution. 
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1. Introduction 

A very fundamental characteristic of galactic disks is the 
flatness, defined by the ratio between exponential scale- 
lengths in the vertical and radial direction h z /h. The 
determination of the thickness of stellar disks is a diffi- 
cult enterprise and has been determined mostly for edge- 
on galaxies, although attempts have been made in other 
cases, too (Ma et al. 1998). Typical values are found to be 
between 0.15 and 1 (van der Kruit & Searle 1981, 1982; 
de Grijs 1998) with a systematic trend for galaxies to be- 
come thinner from SO to Sc if a more complex multicompo- 
nent structure (thick/thin disk) is not taken into account. 
Stellar velocity dispersion measurements show that the 
dispersion tends to grow with total luminosity, and vary 
radially as the square root of surface density, so as to 
support a constant scaleheight (Bottema 1993). Recently 
Kregel et al. (2002) have shown that the flattening of disks 
increases with the amplitude of their maximum rotation, 
and their total HI mass. 

This flatness is an important clue to the formation 
and evolution history of galaxies. Disk stars are thought 
to be formed out of a very thin gaseous disk. The ob- 
served structure then thickens due to scattering by local 
fluctuations of the gravitational field as caused, e.g., by 
giant molecular clouds or spiral structure in the disk. The 
velocity dispersion of the stars increases in both the ra- 
dial and vertical direction by this diffusion through phase 
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space (Wielen 1977, Binney & Lacey 1988). The disk may 
thicken more efficiently through interactions with small 
companions or minor mergers, as suggested by numeri- 
cal simulations (Quinn et al. 1993, Walker et al, 1996, 
Velazquez & White 1999). The thickening could amount to 
a factor of 1.5-2 and can be a constraint on the frequency 
of mergers if disks are observed today to be too thin (e.g. 
Toth & Ostriker 1992). The existence of a thick stellar disk 
in our own Galaxy has been attributed to an old merger 
event (Robin et al 1996, Dalcanton & Bernstein 2002). 
However, other phenomena have to be taken into account 
since the thin disk could also be maintained through gas 
infall with subsequent star formation. Many observational 
facts strongly suggest a high gas infall rate (e.g. Toomre, 
1990, Sancisi et al. 1990, Jiang & Binney 1999). Dynamics 
of galaxy disks, bar and spiral reformation constrain this 
rate such that galaxies may double their mass in less than 
10 Gyrs (Bournaud & Combes 2002, Block et al 2002). 

The influence of tidal interaction on the flatness of 
galaxy disks has been addressed by Reshetnikov et al. 
(1993) and Reshetnikov & Combes (1997) who found 
in interacting galaxies the ratio h z /h to be a factor of 
two higher, due to a disk thicker in absolute values, but 
also shorter in radial dimensions. Schwarzkopf & Dettmar 
(2000, 2001) confirmed this trend on a larger sample, but 
essentially because of a larger absolute thickness of in- 
teracting disks. They also noted that tidally interacting 
galaxies are more frequently warped and that the thick- 
ening was more noticeable in the outer parts. 
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All these studies suggest that the relative flatness of a 
galaxy disk can vary by both effects during its evolution 
since the majority of galaxies have experienced interac- 
tions in the past. It appears that galaxies were dynami- 
cally "hotter" in the past (Abraham et al. 1999), which 
is seen in a smaller fraction of barred galaxies, and this 
could also have some consequences for the flattening. As 
the frequency of interactions/mergers strongly increases 
with rcdshift (e.g. Le Fevre et al 2000), it is interesting 
to trace the evolution of the disk flatness with time by 
studying photometrically distant galaxies with high spa- 
tial resolution. This is now possible in the Hubble Deep 
Fields HDF-N and HDF-S and we present below the re- 
sults of such a study. Section 2 describes our sample and 
Sect. 3 gives the derived general characteristics of disks 
in terms of radial and vertical profiles. The last section 
discusses the main results on the evolution of the flatness 
for stellar disks. 

Throughout the paper, we adopt a flat cosmology with 
O = 1 and H = 70 km s _1 Mpc -1 . (For a model with 
f^m = 1/3, Ha = 2/3 and Cl m + JIa = 1, linear sizes at 
z = 0.7 — 1 will be «25% larger and absolute magnitudes 
are brigthcr by «0^"45.) If not differently specified, all 
magnitudes are expressed in the AB system (Oke 1974). 

2. The sample 

To study the structural parameters of distant edge-on spi- 
ral galaxies, we consider a sample of galaxies presented 
in Paper I (Reshetnikov et al. 2002). This sample includes 
45 galaxies with a diameter larger than l'/3 selected in the 
Hubble deep fields north and south. By applying a V/V max 
completeness test we derived that the sample is statisti- 
cally complete - y/V max =0.50±0.04 (Paper I). From the 
total sample of 45 distant galaxies we rejected objects 
with strong overlaps or interactions (e.g. nl031, SB-WF- 
2353-1914/2340-1898) and objects with unusual morphol- 
ogy (for instance, the "tailed" galaxies nl6, n632, nl027). 
Also, we excluded the galaxies n749 and SB-WF-3329- 
3173 since they seem to be deviating from the edge-on 
orientation. The final sample consists of 34 edge-on galax- 
ies suitable for further analysis. 

The mean isophotal diameter of the galaxies within 
M(^8i4)=26.0 is 2'.'3 (14 kpc at z — 0.9), the mean minor 
axis is 0'.'7. The latter is about 5 times larger than the 
FWHM of the PSF and therefore the galaxies are suffi- 
ciently resolved to study the vertical surface brightness 
distribution. 

The general characteristics of the sample galaxies are 
presented in Table 1. Figure 8 (see Appendix) gives con- 
tour maps and radial major axis profiles of the objects. 

3. General characteristics of edge-on disks 

As judged from to the spectral energy distributions 
(Fernandez-Soto et al. 1999), late-type galaxies (Sbc/Irr) 
predominate in our sample (Table 1). The mean rcdshift 
of the sample and the median of the distribution are the 



Table 1. Edge-on galaxies in the HDF-N and HDF-S 



ID 




b/a 


z 


Type 


15 


24.70 


0.28 


0.45 


Scd 


55 


23.32 


0.20 


0.18 


Irr 


120 


24.27 


0.24 


1.60 


Irr 


153 


24.94 


0.35 


0.73 


Irr 


273 


22.47 


0.45 


0.905 


Ell 


304 


25.60 


0.35 


1.78 


Scd 


450 


24.25 


0.30 


0.70 


Irr 


476 


22.77 


0.39 


0.421 


Irr 


506 


24.47 


0.25 


0.751 


Irr 


534 


23.82 


0.25 


0.95 


Sbc 


671 


23.93 


0.36 


0.681 


Irr 


716 


22.61 


0.30 


0.944 


Sbc 


727 


23.93 


0.33 


0.904 


Scd 


733 


25.04 


0.25 


0.94 


Irr 


774 


21.22 


0.28 


0.485 


Sbc 


805 


25.33 


0.30 


0.66 


Irr 


817 


22.19 


0.40 


0.485 


Sbc 


886 


25.25 


0.30 


0.97 


Irr 


888 


23.70 


0.30 


0.559 


Irr 


898 


25.89 


0.20 


0.94 


Irr 


899 


22.94 


0.20 


0.564 


Scd 


938 


23.03 


0.27 


0.557 


Scd 


979 


23.45 


0.27 


0.517 


Scd 


S0564 


25.27 


0.35 


1.02 


Irr 


S0566 


25.67 


0.30 


2.16 


Scd 


S0578 


23.16 


0.22 


0.98 


Scd 


S0747 


25.31 


0.24 


0.98 


Irr 


S1085 


25.51 


0.29 


1.13 


Irr 


S1404 


23.09 


0.20 


0.50 


Sbc 


S2661 


24.10 


0.30 


0.92 


Irr 


S2691 


25.05 


0.35 


0.46 


Irr 


S3053 


26.07 


0.30 


1.05 


Scd 


S3458 


23.40 


0.40 


0.69 


Scd 


S3685 


26.63 


0.35 


2.23 


SB1 



The columns are: object identification (Paper I); 7si4 magni- 
tude in the AB photometric system according to Fernandez- 
Soto et al. (1999) (aperture magnitude corrected for the effect 
of neighboring objects and the flux lost in pixels outside the 
aperture); apparent axial ratio (from ellipse- fitting of the outer 
contours); redshift (spectroscopic redshifts are with three num- 
bers, photometric redshifts obtained as the mean of redshifts 
presented by Fernandez-Soto et al. 1999 and by Fontana et al. 
2000); best-fit spectral type of galaxy (Fernandez-Soto et al. 
1999). 

same, z « 0.9. The mean rest-frame absolute magnitude 
of the galaxies is relatively low - M(B) « — 18 m (Paper I). 
This is not unexpected since the Hubble deep fields have 
a very small coverage of only ^2 Mpc (comoving) at z ~ 1 
each. At z < 1, the expected number of bright L* galaxies 
in the two fields is only ^20-60 (Ferguson et al. 2000). 
Therefore, assuming random orientation of galaxies, we 
can expect to find a few highly inclined bright spirals at 
z < 1. 

Two of our edge-on galaxies are possible sources of 
sub- mm (SCUBA) radiation (Serjeant et al. 2002). These 
arc the objects 304 and 774 in Table 1. 
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The mean observed axial ratio of the sample objects 
is 0.3 (Paper I, Table 1), a value exceeding the b/a ra- 
tio for local edge-on late-type galaxies. Possible reasons 
are a large intrinsic thickness of distant objects and also 
the influence of the point-spread function (PSF) on the 
apparent thickness of the small (~ 2") HDF galaxies. 

The largest uncertaintiy arises from the range of galaxy 
inclinations in our sample since we do not know the true 
intrinsic flatnesses of the sample objects and we also can- 
not use a standard approach, e.g., on the basis of the 
Hubble formula or from dust lane geometry. We selected 
our sample on the basis of careful visual inspection of the 
HDF frames. Secondary inclination indicators (asymmetry 
of the minor axis profiles, apparent shift of the nucleus lo- 
cation from the symmetry plane) are in the same range as 
in local samples of highly inclined spirals (e.g. de Grijs & 
van der Kruit 1996, Schwarzkopf & Dettmar 2000). Also, 
comparison with results of realistic numerical models for 
inclined spirals (e.g. Byun et al. 1994) indicates that most 
of our sample objects are within 5° of edge-on orienta- 
tion. Such a moderate deviation from edge-on orientation 
does not significantly change the slope of the vertical sur- 
face brightness distribution (van der Kruit & Searle 1981, 
Barteldrees & Dettmar 1994). 

3.1. Optical colors 

Figure [l] compares the observed colors of the 14 sample 
objects with spectroscopic redshifts and of nearly face-on 
spirals (randomly selected from the HDF-N) as a function 
of redshift. Both subsamples demonstrate the same trend 
of colors with z with a somewhat larger scatter for the 
edge-on disks. In the redshift interval z = 0.4 — 1 the 
edge-on disks show a hint of reddening in comparison with 
the face-on spirals: = +0.25 ± 0.18 and Ef,- V = 

+0.23 ± 0.13. These estimates correspond to inclination 
corrections for local spiral galaxies (e.g. de Vaucouleurs et 
al. 1991). 

The predictions of apparent color evolution for two 
models of galaxy formation are shown in Fig. |l| according 
to Westera et al. (2002). The first model is characterized 
by a slowly growing dark halo with a continuous infall of 
gas and dark matter (accretion model, dashed line). The 
second one is a classical collapse model in a dark matter 
halo (collapse model, dotted line) . Obviously, both models 
provide a satisfactory description of general observational 
trends of colors with redshift. 



3.2. Major axis profiles 

To investigate the global structure of the sample objects 
we adopted as a fitting model the simple standard dou- 
ble exponential disk: L(r,z) — LoCxp(— ^ — W), where h 
and h z are the scalelength and scaleheight, respectively. 
This simple model gives a reasonable description of the 
structure for most local disks, even in edge-on orientation 
(e.g. Misiriotis et al. 2000). Comparison with other verti- 
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Fig. 1. Observed color - spectroscopic redshift depen- 
dence for the HDF edge-on galaxies (solid rhombs) and 
for face-on spirals (open circles). The u — b, b—v, and v — i 
corresponding to U 300 - B i50 , B 450 - V 606 , and V 606 - hu 
colors (Williams et al. 1996). The lines show predicted 
color evolution as a function of redshift for two models of 
disk galaxy formation - the accretion (dashed line) and 
the collapse (dotted line) models (Westera et al. 2002). 



cal models for disks (for instance, with "sech 2 (z/z )" or 
"sech(z/zo)" distributions) can be done via h z = z /2 = 



z>/V2. 



3.2.1. Exponential profiles 

Figure 8 presents the major axis photometric profiles of 
all galaxies in the Isu passband. The profiles show a large 
diversity in their shapes. Trying to fit the surface brighness 
distributions by an exponential law, we find that only 19 
or 56%±7% of the disks are reasonably well described by 
exponential-like profiles. 

The parameters for the exponential disks are presented 
in Table 2. Typical errors for the scalelength values are 
w 10%. The errors of scaleheights are presented in Table 
2 and they are not just the formal errors of the least- 
squares fit, rather they were estimated by comparing re- 
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suits from four independent fits (see Sect. 3.3). To check 
the reliability of the exponential model to describe the 
global photometric structure of the studied galaxies, we 
have calculated total luminosities as L, 



cxp 



2nl h 2 b/a, 

where Io=dex(-OA^,( ) ) (Table 2) and b/a is the appar- 
ent axial ratio (Table 1). We found that the mean differ- 
ence of observed and model magnitudes is (Isu ~~ Icxp} = 
+0.31±0.12. This small difference can be attributed to the 
infinite extrapolation for the model magnitudes. 

The mean observational characteristics of the expo- 
nential disks are (Mhu)) = (22.99± f .1 6) mag /D", (h) = 
0"43±0'.'fl, (z) = 1.0±0.5. The observed values of/xo(Isi4 
(uncorrected for any inclination effects) have been con- 
verted to a rest-frame B by applying the cosmological dim- 
ming term and a fc-correction color term for Irr galaxies 
(Lilly et al. 1998, Lilly et al. 1995, Paper I). The resulting 
mean characteristics of the central surface brightness and 
of the scale length (after a small correction for the PSF - 
see Sect. 3.3) are (fi (B)) = (20.7±1.3) mag /D", (h) = 2.3± 
0.7 kpc. For the disks with z < 1 (14 objects), the central 
surface brightness is (21.05±1.0) mag /D". Therefore, dis- 
tant edge-on disks demonstrate somewhat brighter central 
surface brightnesses in comparison to nearby "Freeman" 
disks for which fj, (B) = 21.5 mag /D" in the AB sys- 
tem (Freeman 1970). However, this difference is negli- 
gible within the quoted scatter. Note also that our de- 
rived value of the mean central surface brightness is in 
good agreement with results for local samples of edge- 
on galaxies: (fi (B)) = (21.0 ± 0.6) mag /D" (28 mem- 
bers of interacting systems, Reshetnikov & Combes 1997), 
(Mo(-B)) = (21.0±l.l) mag /D" (45 normal galaxies, deGrijs 
1998). 

As a first order approximation, let us assume that the 
galaxy centers are optically thick, or close to that, and so 
the central surface brightness does not depend strongly on 
the inclination of the disk. This is a quite realistic approxi- 
mation since at redshift z = 0.9 the filter Ign corresponds 
approximately to filter B (8140 A:(l+0.9)=4280 A) and 
the central regions of local spirals are not optically thin 
in the blue spectral region. There are several lines of ar- 
guments - discussed in the following - in favour of such a 
conclusion. 

The only direct and model-independent way to study 
the transparency of the centers of galaxies is the chance 
superposition of galaxies. In the innermost few hundred 
parsecs of NGC 3314A (a sub-L* Sc spiral galaxy) even 
the most transparent regions between dust lanes show an 
extinction of A B w 7 rn (Keel & White 2001). The disk of 
NGC 3314A is inclined by i = 41°, so this galaxy is far 
from edge-on orientation. 

The same conclusion can be derived from statistical 
arguments. For instance, Giovanelli et al. (1995), using a 
sample of more than 1700 Sc galaxies with /-band pho- 
tometry, found only a weak dependence of the central sur- 
face brightnesses of the disks on inclination: /io (I) bright- 
ens by about 0.3-0.4 mag between the face-on and edge- 
on aspects. In the B band the dependence must be even 
much weaker. Also, with regard to distant galaxies, Lilly et 
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Fig. 2. fio(B)-h plane for edge-on galaxies. Typical error 
of the data is shown by a small cross at bottom left. The 
dotted line shows the Freeman law, the solid line is the 
line of constant disk luminosity. The solid cross gives the 
mean characteristics of our edge-on disks and the dashed 
cross is the mean for the sample of 917 local late type 
spirals with M(I) > -22 m from Byun (1992). 



al. (1998) did not find any significant correlation between 
fJ-o(B) and disk axial ratio for spirals with z = 0.2 — 1. 

Numerical simulations of the radiative transfer in 
dusty disks show that even a small optical depth can sig- 
nificantly reduce the expected projection effect on the cen- 
tral surface brightness, especially in the B band (Byun et 
al. 1994). 

Therefore, assuming that the central surface bright- 
ness shows only a weak (or negligible) dependence on the 
galaxy inclination (Afi (B) = (0 - 0.5) mag /D") we can 
estimate an average "face-on" absolute magnitude of the 
disks as (M(B)) = (no(B)) - 51og(ft) - 38.57 = -19™6...- 
1 9™ 1 . (Let us note also that in edge-on orientation the 
presence of dust makes the scalelength appear larger than 
in the face-on case (e.g. Byun et al. 1994) and this can re- 
duce the estimated face-on luminosity.) The observational 
value for edge-on disks is -18™1 or 1™0...1™5 fainter. The 
value of dimming (l^O.-.l^S) is consistent within uncer- 
taintes with the value of global extinction for spirals (e.g. 
de Vaucouleurs et al. 1991, Tully et al. 1998). 

Table 2 presents also the h values in the Veoe pass- 
band. The mean scale length ratio is (h(Isi4) / h(Veoe)) 
= 1.04±0.13. (In the rest-frame this ratio resembles 
h(B)/h(U).) Although we did not find any significant 
color dependence as expected from the modeling by de 
Jong (1996) for these filters, we cannot draw any firm 
conclusions due to the small-number statistics. 

In Fig. H we consider the fJ,o~h plane. The dashed cross 
in the figure shows the mean parameters of 917 Sb-Sd spi- 
ral galaxies from Byun's (1992) thesis (converted to Bab)- 
The sample of Byun is angular diameter- limited (like ours) 
and contains galaxies with angular diameter d > 1.'7. The 
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mean absolute magnitude of this sample of local galax- 
ies is (M(I)) = -20 m 4, or (M (B)) w -18TT (adopting 
B — V=+1.7 for the disks of galaxies, de Jong 1996). 
Therefore, in edge-on orientation the local galaxies will 
show about the same luminosities as our sample galaxies. 
As one can see in the figure, the mean characteristics of 
nearby and distant sub-L* disks are close (within quoted 
uncertaintes) but z ~ 1 galaxies show a small displace- 
ment to higher values of surface brightness and smaller 
sizes. The shift to smaller values of h can be even larger 
due to the above mentioned effect of an apparent increase 
of scalelength in edge-on orientation. 

At least three objects (nl5, n805, and n898) can be 
considered as low surface brightness galaxies. These three 
galaxies show the rest-frame central surface brightness 
H {B) w 22 and are intrinsically faint (M(B) > -17 m ). 
Adopting small or neglible internal absorption for such 
galaxies (e.g. Tully et al. 1998) and converting the obser- 
vational values of i^o(B) to face-on orientation, we obtain 
nl{B) « 23.5 - 24 ma s/D", typical for local low surface 
brightness galaxies (e.g. UGC 7321 - Matthews & Wood 
2001). 
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Fig. 4. Photometric profiles of HDF-N 866. Solid line - 
P.A.=125°, dashed line (the profile shifted down by l m ) 
- P.A.=35°. The P.A. orientation is according to Fig. |^ 
(from top to the left). 



3.2.2. Non-exponential profiles 

15 of 34 disks (44%±5%) show "ridge" -like profiles, with 
several peaks and are often very asymmetric (see, e.g., 
n55, n506, or n886 in Fig. 8). Therefore, roughly half of 
z ~ 1 disks are not exponential. This fraction may be 
even larger due to preselection of the sample (Sect. 2). This 
strongly contrasts with the local Universe and supports a 
late formation of disks (e.g. Mo et al. 1998, see van den 
Bergh et al. 2000 about galaxy morphology evolution with 
redshift). 




Fig. 3. Reproduction of the HDF-N galaxy 866 
(2."3x2."3) in the I su passband. 

To understand the origin of non-exponential brightness 
distributions, we have looked for possible counterparts for 



such objects among face-on galaxies in the deep fields. One 
of the best examples is the HDF-N galaxy 866 (2-153.0 in 
Williams et al. 1996) with the photometric redshift 0.96, 
7gX4 = 23.6 and the spectral type Scd (Fernandez-Soto et 
al. 1999). The observed characteristics of this galaxy are 
close to those for the sample edge-on objects. This galaxy 
exhibits an irregular patchy disk with a low-contrast nu- 
cleus and embedded bright knots (five at least - see Fig.|||) 
and no spiral structure has developed yet. Figure ^ shows 
two photometric profiles of the galaxy which resemble the 
"ridge-like" profiles of many of our edge-on disks men- 
tioned above. The galaxy is very similar in optical mor- 
phology to the peculiar galaxy H36490_1221 discussed in 
van den Bergh et al. (2001). Those authors suggest that 
such objects represent proto-late-type galaxies. Therefore, 
we can speculate that about half of our sample objects 
are young disks still in the formation process which will in 
their future evolution relax to an exponential distribution. 

To check the frequency of non-exponential profiles 
among less distant spirals we selected in two deep fields all 
bright (Isi4 < 22™ 1 ) and non-edge-on galaxies of Sbc-Irr 
types. In total, we selected 17 galaxies with mean red- 
shift (z) « 0.5 and absolute magnitude (M(B)) « -20™ 
Almost one third (5/17) of these objects show strongly 
non-exponential surface brightness distributions (work 
in preparation). The rest of the objects often demon- 
strate significant (~ 0™5) deviations from the exponent. 
Therefore, non-exponential density distributions of galaxy 
disks are indeed very frequent at z > 0. 

3.3. Minor axis profiles 

The typical scalelengths of our objects are » 0'.'3 — 0'.'6 
(Table 2). Therefore, h z values are expected to be < O'.'l 
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Table 2. Characteristics of galaxies with exponential- like surface brightness distribution 



ID 


Mo 


h 1 


V 

Mo 


h v 


hi/h 1 


h 1 


hi ± a hz 






(") 




(") 




(kpc) 


(kp 


c) 


15 


22.94 


0.38 


23.27 


0.46 


0.18 


1.70 


0.31 


0.06 


153 


23.17 


0.29 


23.69 


0.28 


0.21 


1.50 


0.32 


0.10 


304 


23.70 


0.32 


24.05 


0.35 


0.27 


1.73 


0.47 


0.25 


450 


23.34 


0.48 


23.79 


0.48 


0.33 


2.62 


0.86 


0.34 


476 


21.67 


0.40 


22.22 


0.46 


0.27 


1.74 


0.47 


0.09 


671 


22.54 


0.39 


23.25 


0.37 


0.23 


2.03 


0.47 


0.16 


727 


22.28 


0.36 


22.50 


0.34 


0.18 


1.98 


0.36 


0.06 


733 


23.76 


0.59 


24.54 


0.70 


0.17 


3.32 


0.56 


0.11 


774 


20.42 


0.31 


21.57 


0.32 


0.19 


1.40 


0.27 


0.03 


805 


23.56 
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0.21 
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0.41 
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0.03 
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0.34 


0.14 
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24.00 


0.28 
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0.24 


0.22 


1.43 


0.31 


0.11 



The columns are: galaxy identification; central surface brightness and exponential scale length in arcsec (isi4 passband); the 
same for the Veoe filter; scale height to scale length ratio corrected for the PSF; final values of scale length and scale height in 
kpc. 



or below the PSF width of FWHM(PSF)=0'/14. It is im- 
portant to note that h z values are characteristics of the 
surface brightness gradients and are not linear sizes of any 
structures. To determine true scaleheights, we must study 
the effects of the HDF's PSF on the surface brightness 
distribution of edge-on galaxies. 

We have created a grid of double exponential disk 
models with various input parameters in the region h = 
0"3 - 0"6 and h z /h from 0.10 to 0.35. Then we prepared 
an empirical PSF by averaging several stellar images from 
the /8i4-band frames of the HDF and convolved the model 
disks with the PSF. Comparing input and output disk 
characteristics, we determined empirical corrections for 
the observed parameters. The output parameters were es- 
timated from the analysis of major and minor axes cuts 
of the model images. Finally, we found a strong influence 
of the PSF on the vertical scale of a disk. Also, we noted 
a small (within 10%) overestimation of the scale length 
values due to the PSF. 

To present the influence of the PSF more clearly, we 
show in Fig. [s] the dependence of vertical stretching due 
to the PSF on the intrinsic flattening (h z /h) of the model 
(the factor S gives the ratio of observed and true values of 
flattening). It is evident that intrinsically small (h ss 0'.'3) 
and thin (h z /h sa 0.1) galaxies become more than 2 times 
thicker after application of the PSF. Therefore, we cannot 
expect to find very thin galaxies among small and distant 
HDF objects. 

We excluded the central regions of the galaxies from 
the study in order to avoid any possible contribution of 
the bulges. (Although we do not see clear signs of bulges 
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Fig. 5. Dependence of the "stretching" factor S due to 
the PSF on the true scaleheight to scalelength ratio h z /h. 
The lines marked as 0.3, 0.4, 0.5 correspond to models 
with h=Qf.% 0'.'4, 0'.'5. 

present in most galaxies - see Fig. 8). In order to explore 
the vertical structure of 19 exponential objects we exam- 
ined one-dimensional cuts extracted at ±(0.5-l)/i in the 
galactocentric distance, avoiding the faintest outer regions 
of the objects. To minimize the influence of dust, we ex- 
cluded the 1-2 central pixels of the cuts. In the analysis, 
we distinguished between the different sides of the galaxy 
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planes and obtained the mean value of h z for each galaxy 
from 4 independent estimates. The mean observational 
(uncorrected for the PSF influence) scale height for the 
sample objects is 0'.'13, the mean dispersion a hz = 0'.'03. 
This value of dispersion (23%) characterises the combined 
effect of the intrinsic variations of vertical scales and of 
errors in scale estimates. 

The results of the analysis (corrected for the PSF) are 
given in the last three columns of Table 2. The mean rela- 
tive thickness of the galaxies reduces after the PSF correc- 
tion from (h z /h) obs = 0.31 ±0.09 to (h z /h) = 0.21 ±0.08. 

3.3.1. Vertical scale 

The mean value of the vertical scale for distant edge-on 
galaxies is (/i z )=0.46±0.21 kpc or (z ) « 0.9 kpc. The 
obtained mean value is unbiased since we work with a 
statistically complete sample of galaxies with angular di- 
ameter larger than 1'.'3 (Sect. 2, Paper I). These values are 
typical for normal spiral disks at z « (e.g. Reshetnikov 
& Combes 1997, Schwarzkopf & Dettmar 2000). However, 
the direct comparison of local and distant samples must 
be carried out with caution due to very different selec- 
tion effects. For instance, distant galaxies are sub-L* and 
relatively compact while the local samples usually con- 
tain brighter and larger disk galaxies (see Table 2 in 
Reshetnikov & Combes 1997). 

In order to avoid (at least partially) the selection ef- 
fects, we decided to compare our sample to local galax- 
ies of the same size. Figure ^|a shows the dependence of 
scaleheight on scalelength for 38 edge-on spiral galaxies 
(T > 0) in the I passband according to de Grijs (1998). 
He adopted an exponential model for the observational 
brightness distribution and considered one-dimensional 
photometric profiles so we can directly compare our re- 
sults. The mean regression for local galaxies predicts that 
a galaxy with h = 2.3 kpc must have h z = 0.3 kpc. 
Therefore distant galaxies are thicker by a factor 1.5 on 
average. (In the non-zero A cosmology with VL m = 1/3, 
£1a = 2/3 the relative thickening stays the same.) 

Figure ||b presents the results of a two-dimensional 
decomposition in the / passband of 34 edge-on spirals 
from the de Grijs (1998) sample recently presented by 
Kregel et al. (2002). The fitting model assumes a trans- 
parent disk (this is very questionable for our sample ob- 
jects) and includes the line-of-sight integration of the lu- 
minosity distribution. Comparison of two methods shows 
that a one-dimensional method can overestimate the scale- 
lengths by about 20% (Figure 3 in Kregel et al. 2002). 
However, a direct comparison of h values shows that h(de 
Grijs)//i(Kregel et al.)=0.93±0.20. A comparison of scale- 
height values obtained by the two methods leads to a ratio 
h z (de Grijs) /h z (Kregel et al.)=1.13±0.22. Strictly speak- 
ing, both methods give a statistically insignificant differ- 
ence. The Kregel et al. (2002) analysis results in a ~1.4 
thickening of distant disks. 



In Figure gc we compare the Schwarzkopf & Dettmar 
(2000) sample of non-interacting galaxies (the data in the 
Johnson R or Thuan & Gunn r filters) . The fitting model 
avoids the central dust lane in the line-of-sight integration. 
In this case the thickening factor increases to 1.9. 

The last sample (Figure ^d) includes 153 edge-on spi- 
ral galaxy with published photometric parameters in the 
K s band (Bizyaev & Mitronova 2002). The fitting model 
assumes a transparent disk. This work gives about a 1.8- 
fold thickening of distant galaxies. 

3.3.2. h/h z ratio 

The relative thickness - h/h z - gives more reliable in- 
formation about the possible thickening. For our sam- 
ple we have (/i//i z )=5.4±2.0, or, excluding one object 
for being most likely not exactly an edge-on galaxy 
(n817), 5.55±1.9. Scaling this in terms of zq, we obtain 
{h/ Zo)=2-7-2.8, with a median value of 2.63. The disks of 
local late-type galaxies are distributed over a wider range 
with the mean h/zo of ss 5 (see Table 2 in Reshetnikov & 
Combes 1997). Moreover, disks of more late-type and gas- 
rich galaxies are thinner (Reshetnikov & Combes 1997, de 
Grijs 1998, Schwarzkopf & Dettmar 2000). For Sc/Sd disks 
Schwarzkopf & Dettmar (2000) give h/z sa 8±2. De Grijs 
(1998) and Kregel et al. (2002) give somewhat smaller ra- 
tios of 3.7 and 4.3, respectively. The Bizyaev & Mitronova 
(2002) sample shows the mean ratio of h/zo — 4.8. 

In the previous discussion we neglected the possible 
dependence of scalelength on the passband. It is known 
(e.g. Table 7 in de Grijs 1998) that the h value increases 
to the blue spectral region. The ratio of scalelength val- 
ues in the B and K passbands can reach « 1.5, and 
h(B)/h{I) w 1.3 (de Grijs 1998). This effect, however, 
would only strengthen the thickness difference in the sam- 
ples of local and z ~ 1 galaxies. 

Therefore, z ~ 1 disks are definitely thicker in compar- 
ison to local galaxies. The factor for the vertical thickening 
(>1.5) is approximately the same as for z m interact- 
ing spiral galaxies (~1.5-2: Reshetnikov & Combes 1997, 
Schwarzkopf & Dettmar 2000). 

The dependence of h/h z values on the rest-frame cen- 
tral surface brightness is shown in Fig. f^. Thin disks are 
present for fJ,o(B) > 21 mag /D", while among bright ob- 
jects (/Uo(-B) < 21 mag /D") we observe only thick disks with 
h/zo 2.2. Taking into account that the -B-band surface 
brightness characterises the current rate of star formation, 
one can speculate that enhanced star-formation and en- 
hanced thickness are the results of ongoing interactions 
and mergings. The observed morphology of galaxies does 
not contradict such a suggestion (Fig. 8). 

3.4. Vertical scale for non-exponential galaxies 

We examined the vertical profiles for 15 non-exponential 
galaxies at distances where the surface brightness drops by 
l m relative to the nucleus. The mean observational value 
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of the scaleheight obtained ((h z ) = 07125) is the same as 
for exponential objects (Sect. 3. 3) with a mean dispersion 
of individul values of 0'.'03. Taking into account that the 
redshift and angular diameter distributions in both sub- 
samples are approximately the same, one can conclude 
that non-exponential disks must have approximately the 
same corrected value of the scale height (h z ^0.5 kpc). 

4. Discussion and conclusions 

Our sample of distant edge-on galaxies shows an unex- 
pectedly large fraction of non-exponential disks. From this 
observation we infer that many sub-L* disks are still in the 
formation process at z ~ 1, or the degree of tidal pertur- 
bations is higher, which has already been observed (e.g. 
Le Fevre et al. 2000). One of the possible mechanisms to 
redistribute matter and to develop the exponential profile 
is the action of a long-lived bar (e.g. Combes et al. 1990). 
The young galaxies of our sample may not have experi- 
enced enough bar torques for this process, either because 
of very transient bars in very unstable disks, or because 



the effects of bars are continuously perturbed by mergers 
and accretion. The consequence is that a galaxy will spend 
much less of its time in a barred state. This conclusion is 
in good agreement with the general depletion of barred 
spirals at z > 0.7 (Abraham et al. 1999, van den Bergh 
2002). 

We have taken into account the effect of the PSF to 
demonstrate that the distant disks examined in this work 
are thicker by a factor of 1.5 if compared to local disk sam- 
ple. This is a preliminary result based on a small number 
of sub-L* galaxies, and should be confirmed by larger sta- 
tistical samples. 

Since the star formation rate increases strongly with 
redshift (e.g. Madau et al. 1996, Flores et al. 1999, 
Thomson et al. 2001) one could speculate that the fraction 
of young stars formed in the midplanc of disks is expected 
to be much higher in the past, making the young disk ap- 
pear much thinner. Previous work on local samples has 
already concluded that tidal interactions between galax- 
ies have a certain thickening effect of the order of 1.5-2 
(Reshetnikov & Combes 1997, Schwarzkopf & Dettmar, 
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Fig. 7. Dependence of the relative thickness of distant 
galaxies on their rest-frame central surface brightness. 

2000). Since the frequency of galaxy interactions/mergers 
increases strongly with redshift, the present result fits into 
this picture. If one assumes that all disks have undergone 
such interactions, either additional cooling processes for 
the hot disks or a substantial growth predominantly in 
the radial direction are required to explain the flatness of 
disks in the local Universe. 

However, other parameters should be taken into ac- 
count in the interpretation. The distant disks are on ab- 
solute scales smaller than the local sample disks. It is 
likely that the morphological types are later (less evolved) 
and the masses of the galaxies smaller on average, as al- 
ready found in previous studies based on HST surveys 
(e.g. Abraham et al. 1996, Glazebrook et al. 1998). The 
lower frequency of bars, associated with the smaller bulge- 
to-disk ratio, suggests that disks are more unstable and 
dynamically hotter (Abraham et al. 1999), which could 
explain the present thickening effect. In this context it is 
also noteworthy to mention that the galaxies of our sam- 
ple contain only a few bulges as can be judged from the 
major axis profiles in Fig. 8. This can either be interpreted 
as evidence for a later formation of bulges by secular pro- 
cesses (e.g. Combes 2001) or that our sample is selected 
such that it represents early phases of late morphological 
types only. 
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5. Appendix 

We present here contour maps and major axis photometric 
profiles for the sample of edge-on galaxies in the HDF (Table 
!)• 




Fig. 8. Contour plots (left) and major axis photometric profiles (right) of edge-on galaxies in the Isu passband. The 
faintest contour corresponds to 2a limiting isophote of the HDF or /i(isi4) ~ 25.75, isophotes step is 0. m 5. Dashed 
lines - symmetric exponential model for the disk. 




Fig. 8. (cont.) 
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